Glassiness in Uniformly Frustrated Systems 
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We review several models of glassy systems where the randomness is self generated, i.e. already 
an infinitesimal amount of disorder is sufficient to cause a transition to a non-ergodic, glassy state. 
We discuss the application of the replica formalism developed for the spin glass systems to study 
the glass transition in uniformly frustrated many-body systems. Here a localization in configuration 
space emerges leading to an entropy crisis of the system. Using a combination of density functional 
theory and Landau theory of the glassy state, we first analyze the mean field glass transition within 
the saddle point approximation. We go beyond the saddle point approximation by considering the 
energy fluctuations around the saddle point and evaluate the barrier height distribution. 



PACS numbers: 75.10.Nr, 05.70.Ln 

I. INTRODUCTION 

The glass transition and slow glassy dynamics are 
phenomena most widely studied in the context of su- 
percooled liquids and polymer melts [TJ [2]. Important 
progress has been made in understanding the complexity 
of these phenomena by combining dynamical approaches 
with the concept of an underlying energy landscape, as 
is discussed in detail in other chapters of this volume. 
For example, on the level of the mean field theory of 
glasses, it has been demonstrated that the dynamical, 
ideal mode coupling theory [3j and energy landscape 
based replica mean field theories describe the same un- 
derlying physics, yet from rather different perspectives 
[H [5] . Novel replica approaches have been developed to 
characterize the emergence of a metastable amorphous 
solid [16lll9| . In addition, droplet arguments have been 
proposed to include important physics beyond the pure 
mean field description and led to the formulation of the 
random first order transition theory of glasses 20J . The 
latter is particularly important if one wants to make spe- 
cific predictions for experiments that require one to go 
beyond the mean field limit [2"TH2"5] . This random first 
order transition theory of glasses with an underlying en- 
tropy crisis [26 offers a general and quantitative descrip- 



tion of structural glasses. In addition, it defines a uni- 
versality class for complex many body systems that is of 
importance for a much larger class of materials. 
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FIG. 1: Free energy of a system with entropy crisis as a func- 
tion of the generalized coordinate q. For a specific case of a 
structural glass q plays a role of density of the local atomic 
configurations. Glassy behavior is characterized by the emer- 
gence of large number of metastable states separated by high 
energy barriers. 

Glassy dynamics usually occurs when a system is su- 
percooled below a first order transition to an ordered 



state. In supercooled liquids this ordered state is the 
crystalline solid. If one applies the same ideas to sys- 
tems with charge or spin collective degrees of freedom in 
a correlated material, the ordered state corresponds to an 
electron crystal or an ordered magnetic state. The equi- 
librium state is therefore ordered. Supercooling the liq- 
uid state is possible once the nucleation barriers are high. 
Thus, in what follows we ignore this ordered state and 
consider time scales that are short compared to the nucle- 
ation time t nuc i. One can now ask the question whether 
and for how long one can still consider the supercooled 
liquid in local equilibrium. Here, by local equilibrium 
we mean that the system is ergodic and all configura- 
tions are being explored with a probability given by the 
Boltzmann distribution (only excluding the sharp peak 
that corresponds to the ordered state) . One expects that 
this supercooled state will fall out of equilibrium over 
very long time scales, i.e. behave nonergodically, once 
the barriers between distinct configurations of compara- 
ble energy become very large, Fig. [T] If this happens, the 
system will explore only a subset of the available phase 
space and will eventually freeze. Suppose there are Af ms 
of such metastable configurations, that are separated by 
high barriers, Fig. [2] A calorimetric measurement will 
then yield an entropy that is reduced by 

S c ~ k B logJV ms (1) 

compared to the ergodic situation, where S c is referred to 
as the configurational entropy. If Af ms is exponentially 
large in the system size, S c becomes extensive and the 
dynamically frozen state strictly speaking does not obey 
the third law of thermodynamics. For liquids it appears 
that S c on infinite time scales would extrapolate to zero 
at a temperature Tk > 0. Extrapolating further to T < 
Tk would give S c less then zero. This is the entropy crisis 
of Kauzmann [16]. Such a crisis is avoided by an ideal 
glass transition in the models highlighted here. 

A crucial question to explore is under what circum- 
stances the configurational entropy is a well defined and 
meaningful concept. Metastable states become sharp 
conceptually only within a mean field theory description, 
where barriers can be infinitely high and a system might 
get trapped in a local minimum of phase space forever. 
Within mean field theory, the number of states that are 
separated by high barriers proliferates at the mode cou- 
pling temperature Ta- The sudden onset of S c (T) at 
T = Ta should not be misunderstood to mean that new 
configurations emerge at that temperature. These con- 
figurations were present already in the liquid state above 



Ta- What happens at Ta is that the barriers grow and 
that the system's dynamics becomes localized in configu- 
ration space. It is this localization in configuration space 
that is the key element of the random first order theory 
(not the ultimate glass transition). Within mean field 
theory the system freezes right at T4, while activated 
events, that are a key element of the theory, allow for a 
slowed down dynamics even below Ta, reflecting the fact 
that those barriers are large but finite. 

One might criticize this approach since going beyond 
mean field theory does not allow for a sharp thermody- 
namic definition of S c . However, there are well known 
examples that demonstrate the opposite: understand- 
ing overheated or supercooled phases close to a first or- 
der transition also starts from the mean field concept of 
a local minimum in the energy landscape. Once sup- 
plemented with an appropriate droplet theory of the 
metastable state it yields a description of the nucleation 
dynamics at first order transitions [27] ■ The logic used 
in the random first order theory of glasses is very simi- 
lar. One first develops a mean field theory and then sup- 
plements it by an appropriate droplet calculation that 
clearly goes beyond strict mean field theory. 
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FIG. 2: Emergence of an exponentially large number of 
metastable configurations and the possibility for a system to 
realize these configurations gives rise to the configurational 
entropy. The latter serves as the driving force for the struc- 
tural transitions in the glassy phase. 

There are numerous indications that the underlying 
principles that govern glassy behavior in supercooled liq- 
uids apply to other many body systems as well. Gener- 
ally, glassy systems are characterized by slow relaxations 
and a broad spectrum of excitations, a behavior found 
in a number of strongly correlated hard condensed mat- 
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ter physics systems. A particularly interesting situa- 
tion occurs when interactions, that change very differ- 
ently with distance, compete with one another. In such 
uniformly frustrated systems an intermediate length scale 
can emerge that leads to new spatial structures and inho- 
mogeneities. A classical example is a ferromagnet where 
finite range exchange interactions compete with the long 
range dipole-dipole interactions and leads to the forma- 
tion of magnetic domains. Similar behavior is in fact 
abundant in a number of correlated electron systems. 
Examples are stripe formation in doped Mott insulators 
[28H3T)] , defect formation in two-dimensional electron sys- 
tems [10H32] , bubbles of electronic states of high Landau 
levels in quantum Hall systems [43J 04], magnetic do- 
mains in magnetic multilayer compounds 45J , and meso- 
scopic structures formed in self-assembly systems |46j . 
These systems typically exhibit a multi-time-scale dy- 
namics similar to the relaxation found in glasses. This 
observation suggests that glassy behavior and large relax- 
ation times are caused by the competition of interactions 
with different characteristic length scales [12] I47H56] , and 
not primarily due to the presence of strong disorder in 
the system. As we will see, this does not mean that dis- 
order is irrelevant for uniformly frustrated systems. In 
fact, the opposite is true. We find that even the smallest 
amounts of disorder and imperfections drive such sys- 
tem into a non-ergodic regime that otherwise emerges 
only in the presence of very strong disorder. In the lan- 
guage of a renormalization group flow this means that the 
strength of disorder, g, flows to larger and larger values 
and reaches a strong-disorder fixed point <?* , even if the 
physical, bare disorder strength is very small. Glassiness 
arises spontaneously even for infinitesimal extrinsic dis- 
order, leading to self- generated randomness. Fig. [3] In 
uniformly frustrated systems self generated randomness 
and the emergence of an entropy crisis were first discussed 
in Ref. [S3]. For a beautiful example of self generated 
randomness in the context of coupled Josephson junction 
arrays, see Ref. [57] . 

While frustration through competing interactions is 
easiest to analyze, the notion that uniform geometric 
frustration is at the heart of the structural glass tran- 
sition was already discussed and analyzed in Refs.jBHS] 
with the underlying view that one can capture frustration 
through an underlying ideal structure in curved space as 
a reference state. Recently this general idea was ana- 
lyzed in the specific context of glass formation in a hy- 
perbolic space[5]. Modeling frustrated interactions by a 
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FIG. 3: Phenomenon of self-generated randomness appears in 
a variety of physical systems ranging from strongly correlated 
two-dimensional electrons to polymer melts and mesoscopic 
self-assembly structures. 



competing interactions on different length scales was sug- 
gested in Ref.pUj where the concept of avoided critical- 
ity for uniformly frustrated systems was introduced, see 
also Refs. [HUE] f° r further details. Evidence for fragile 
glass-forming behavior in the relaxation of Coulomb frus- 
trated three-dimensional systems by means of a Monte 
Carlo simulation was given in Ref.[13j. Furthermore, in 
Ref.[H] the close relation between the replica approach 
presented here and the results obtained within a mode 
coupling analysis of a uniformly frustrated system, was 
demonstrated explicitly. For a recent review on the sub- 
ject, see Ref. [TSJ 



In this chapter we discuss in some detail model sys- 
tems that display self generated randomness. We do so 
by applying the replica formalism originally developed 
to describe disordered spin glasses [TB] and structural 
glasses [171 HI]- We demonstrate that there are many 
body systems that undergo a mean field glass transition 
for arbitrarily weak disorder [Sj?-56 ■ We then develop 
a Landau theory of the glassy state that reproduces the 
key physical behavior of these uniformly frustrated sys- 
tems [SSI EH] • The appeal of the Landau theory is that 
it easily permits for a generalization beyond the mean 
field limit, where we include instanton events to describe 
dynamical heterogeneity in glassy systems. 
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II. UNIFORMLY FRUSTRATED SYSTEMS: A 
MODEL HAMILTONIAN APPROACH TO GLASS 
FORMATION 

We summarize the main featured of a simple model 
that exhibits glassy behavior due to infinitesimal amount 
of disorder [53l456j . The model exhibits competition of 
interactions on different length scales and there are no ex- 
plicitly quenched degrees of freedom. Consider </?(x): the 
spatially varying amplitude that characterizes the collec- 
tive degrees of freedom of a many body system. Depend- 
ing on the problem under consideration (p (x) corresponds 
to the magnetization m (x) (in case of magnetic domain 
formations) or a density fluctuation relative to the mean 
density, p (x) = po+ip (x). Here, p (x) stands for the elec- 
tron density in case of a doped Mott insulator [T^l [53] 
or for the density of amphiphilic molecules in case of a 
microemulsion or a copolymer [60H65j . 

We consider the effective Hamiltonian 



T~L — T~Lq + Hi, 



(2) 



where Hq contains all terms of second order in <^>(x) 
Specifically, we consider 



(tW 2 (x) + [Vy>(x)] 



+ -J d d x J d d xV(x)y(|x-x'|)^(x'). (3) 

For the non-linear interaction term we assume the local 
form 



where we use 



H in t = J d d xW(p (x)) 



(4) 



(5) 



In Eq. ([3]), ro < favors local order of y(x). Thus, 
in the absence of the long distance coupling V (x — x') 
the system is expected to be homogeneously ordered. 
The gradient term favors homogeneous configurations of 
<p (x) . It is written by taking the the continuum limit of 
an underlying short range interaction. In the case where 
ip (x) refers to the deviation from the mean density, the 
average (p (x)) vanishes by construction. Still the ex- 
pected low temperature minimum of the free energy cor- 
responds to configurations where regions of macroscopic 
size have constant field values (tp (x) ~ const.), sepa- 
rated by domain walls of minimal area, due to the pres- 
ence of the gradient term. This is in fact a toy model 



for macroscopic phase separation. The corresponding or- 
dering temperature can be estimated within mean field 
theory. For d = 3 it follows, for example, T® = 2 „ 4 ^ , 
with momentum cutoff A of the order of an inverse lattice 
constant. 

The situation changes significantly once we include an 
additional long range interaction 



V{x) = l { 



-(d+2) (loY 



(6) 



Here Iq is the new typical length scale that results from 
the competition of the gradient term and the long range 
interaction with < r < d. The long range interaction 
vanishes in the limit Iq — > 00. The exponent r determines 
the rate of decay of the interaction. The effect of this long 
range interaction can easily be seen in momentum space. 
The Fourier transform of V (x) is given as 



V (k) = C'd- 



r 2 
'0 



sd-r ■ 



(7) 



(l k) 

with dimcnsionless coefficient Cd = 
2 d-T n d/2 T (<*=Z) /r (t/2). For t < d, homogeneous field 
configurations with k — > are energetically very costly 
and spatial pattern with finite wave number kg ~ 2k /Iq 
emerge as consequence of the competition between the 
short range interaction and long range forces. For r = d 
, V (fc) oc log (klo) and homogeneous configurations are 
still suppressed, albeit only logarithmically. This is 
relevant for the dipole-dipole interactions in d = 3. 

One can also consider situations with a screened inter- 
action, such as 



V scr (x) = l (d+2) f ^ exp (-x/a scr ) 



(8) 



As shown in Ref. [63 , glass formation is virtually un- 
changed compared to the unscreened interaction once 
Iq <C ft SC r, while it is suppressed in the opposite limit. 
The theory of Refs. [Ml [S3] was also applied to study 
the glassy phase in a model for charged colloids [55] [57] . 

At the Hartree level, the correlation function that re- 
sults from these two competing interactions is 

1 



G H (k) 



k 2 +r> + V (k) ' 



(9) 



with r' = r + u 4 T J Gh (k). The intermediate 
length scale that results from the competition between 
short and long range interactions leads to a peak in the 
correlation function at 



ko = Cdl 



d l Q 



(10) 
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where Ca 



(d-r) 



We can approximate 



Gn(k) close to fco as: 
G H (k) 



A 



(ii) 



where Aq = (2 + d — r) . The liquid state correlation 
length is given by = A r' + Its temperature de- 
pendence must be determined from a conventional equi- 



librium calculation. The form Eq. (Ill holds for a large 
class of models with competing interactions. For 113 = 
Brazovskii showed that such a model undergoes a fluctua- 
tion induced first order transition to a state with lamellar, 
or stripe order [68]. In general, this model has been dis- 
cussed in the context of complex crystallization [551 [75] • 
For example, Alexander and McTague [70] argued that 
crystallization of body centered cubic crystals is preferred 
if the first-order character is not too pronounced. As 
shown in Refs. [7TJ [72] the preference for bcc order is 
rather for metastable states that form near the spinodal. 
Below we will ignore ordered crystalline states and as- 
sume that their nucleation kinetics is sufficiently slow to 
supercool the disordered state. 

Finally, we mention that our model, Eqs. ([2][4]), can 
also be used to describe interacting liquids consisting of 
A-components [55]. In this case ip — (ipi, • • • , ipjsr) is a 
vector whose components refer to the density deviations 
of the different components of the liquid, with mean den- 
sity po,z- Now, the Gaussian part 



Ho =~ I d d xd d x' 



^i^(x) c w (|x-x'|) +£«<p j <fi> (x') 



(12) 



is determined by the direct correlation function c«< (x) 
of the fluid [72] . Within the density functional approach 
pioneered by Ramakrishnan and Youssoff [75] the non- 
linear part of the Hamiltonian is determined by the ideal 
gas free energy: 



(am + w) 



log 



Po.i + <pi 
Po,i 



2po,i 
(13) 

In case of the density functional theory, Eq. ^ corre- 
sponds to the resummed virial expansion of the liquid. 

In summary, models of the form Eqs. ([2j[4j) capture 
the physics of a large class of systems where competing 
interactions on different length scales lead to complex 
spatial pattern and potentially to slow glassy dynamics. 
In the next section we discuss in some detail the methods 



that will be used to describe and analyze glassy systems. 
Those methods will then be applied to the model of Eqs. 
S2M further below. 



ENTROPY CRISIS AND MEAN FIELD 
FORMALISM 



III. 



We first outline the main idea of the generalized replica 
mean field approach for systems with an entropy cri- 
sis [T5] [T7] . We start from a classical field theory with 
Hamiltonian, H [ip], and field variable, <£>(x), that yields 
the partition function 



D<pexp(-PH[<P])- 



(14) 



The thermodynamic free energy is then given as F = 
—T log Z. Suppose we know the thermodynamic density 
of states 



/'(/)= / DsHi-^m 



(15) 



where V is the system volume. It follows with s c (f) = 
V logp(/) that 

Z = y /mM df exp [-V (/?/ - s c (f, T))] . (16) 

For large V the integral is evaluated at the saddle point, 
leading to the free energy density 

F = -T log Z « V mm [/ - Ta c (f)] . (17) 

In the case when the number of distinct configurations 
with given energy / is exponentially large, s c (/) will be 
finite as V — > 00. Due to the configurational entropy 
density, s c , the free energy density then differs from its 
minimum value, f min , even at the saddle point level. In 
mean field theory s c (T) is finite below a temperature Ta 
that coincides with the dynamic mode coupling tempera- 
ture [J]. At a lower temperature, Tk < Ta, the entropy 
density vanishes like 



s c {T)^T-T K + o({T~T K )- 



(18) 



and the system undergoes a transition into a statically 
frozen state. 

In order to have an explicit expression for s c , it is con- 
venient to introduce the partition function [171 118] 

Z(m) = dfexp [-V ((3mf - s c (f, T))] , (19) 
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leading to the free energy density 



T 



/(«*) = -log Z(m) « min [/ - Ts c (f, T)/m] , (20) 

mv f 



such that / 



dmf(m) 

dm 



gives the most probable value 



of / which might also be written as 

f=zj df f exp [ ~ V (/3/ _ Sc(/ ' T))] • (21) 

Furthermore, the configurational entropy density follows 
from 



1 df(m) 



T dm 



(22) 



Next we will describe the ways of calculating the config- 
urational entropy density s c . 

A. Replica formalism 

A systematic approach to calculate F(m) — V f (m) 
was developed in Ref. [TB] . The procedure of Ref. [TB] 
also reproduces the correct results for the configurational 
entropy in several random spin systems. Furthermore, it 
allows for a rather transparent motivation for the intro- 



duction of the variable m of Eq. (20). In what follows, 



we summarize the main idea and some technical steps of 
this formalism. 

We consider the partition function in the presence of 
a bias configuration ip (x) : 



Z[ip] 



Dipe 



(23) 



Here, J Dip... corresponds to the statistical sum over all 
density configurations of the system. At the end of our 
calculation we will take the limit g — > 0, however only 
after we take the thermodynamic limit. For any finite g, 
the configuration ip (x) enters the problem as a quenched 
degree of freedom, analog to random fields in disordered 
spin systems. We note that in distinction to the usual 
random field problem, the probability distribution func- 
tion of metastable configurations P [ip] is not expressed 
in terms of uncorrelated random numbers, but instead is 
determined by the partition function Z [ip] (see below) . 
The free energy for a given bias configuration is 



F [ip] = -T log Z[ip] 



(24) 



Physically F [ip] can be interpreted as the free energy for 
a metastable amorphous field configuration ip (x). Z [ip] 



in Eq. (23) is obviously dominated by configurations 
if (x) ~ ip (x) that correspond to local minima of H [<p] . 
In the replica formalism, no specific amorphous config- 
uration need be specified in order to perform the calcu- 
lation. Rather, the assumption is made that the proba- 
bility distribution for metastable configurations is deter- 
mined by F [ip] according 



P[iP] cxexp(-/3 cff FM) 



(25) 



and is characterized by the effective temperature 
^B^eff = Pcfi — ^bT. This allows one to determine 
the mean free energy 

F = J DipF[i>]P[ip] (26) 

and the corresponding mean configurational entropy 



S c 



DiPP [ip] log P [ip] . 



(27) 



It is physically appealing then to introduce the free 
energy difference, SF, via 



F = F-6F, 



(28) 



where SF gives the amount of energy lost if the system 
is trapped into locally stable states and hence not able 
to explore the entire phase space of the ideal thermody- 
namic equilibrium. If the limit g — ¥ behaves perturba- 
tively, SF = 0. This indicates that the number of locally 
stable configurations stays finite in the thermodynamic 
limit, or at least grows less rapid than exponential with 
V. In this case all states are kinetically accessible. On 
the other hand, if the limit g — > does not behave per- 
turbatively, it means that the number of locally stable 
states, JV ms , is exponentially large in V. This allows 
us to identify the difference between the equilibrium and 
typical free energy as an entropy: 



SF = TS r 



(29) 



The configurational entropy, S c = logA/" ms , is extensive if 
there are exponentially many metastable states. Within 
mean field theory, where the barriers are infinitely high, 
the emergence of S c renders the system incapable of ex- 
ploring the entire phase space. S c is then the amount 
of entropy which the system that freezes it into a glassy 
state appears to lose due to its nonequilibrium-dynamics. 

The approach of Ref. [TH] was successfully used to 
develop a mean field theory for glass formation in super- 
cooled liquids [HI [18] yielding results in detailed agree- 
ment with earlier, non-replica approaches |761 177j . The 



G 



mean values F and S c can be determined from a repli- 
cated partition function 



Z[m} = J D i>Z m [if,] 



via F = 



S-mF (m) and S r 

am v / L 



T dm 



F (m) with 



F(m) 



T 



log Z (to) 



(30) 



(31) 



and replica index to = . Inserting Z [ijj] of Eq. ( 23 ) 



into Eq. (30) and integrating over tp, one gets 

-^5Dr=i H[v a }+g T.T,b=i f d d xip a (x)ip b (x) 



Z TO = 



D 



ipe 



(32) 

which has a structure similar to a conventional equilib- 
rium partition function. The ergodicity breaking field 
tp causes a coupling between replicas which might spon- 
taneously lead to order in replica space even as g — » 0. 
This spontaneous coupling between different copies of the 
system is then associated with a finite S c and thus glassi- 
ness. 

With T c ff = T/to, expressions for F and S c follow from 
the replicated free energy F (T e g) via: 

n2 d(F(T cS )/T cS ) 



F 



-T, 



Sr = - 



cS dT eS 
dF(T cS ) 
dT cS • 



(33) 



These results are in analogy to the usual thermodynamic 
relations between free energy (F — > F(T c ff)), internal 
energy (U — > F), entropy (S — > S c ) and temperature 
(T —> T c ff), see also Ref. [78]. If the liquid gets frozen 
in one of the many metastable states, the system cannot 
anymore realize its configurational entropy, i.e. the mean 
free energy of frozen states is F, higher by TS C if com- 
pared to the equilibrium free energy of the liquid, Fig. [4] 

The physically intuitive analogy between effective tem- 
perature, mean energy, F and configurational entropy 



to thermodynamic relations, Eq. (33), suggest that one 



should analyze the corresponding configurational heat ca- 
pacity 



C r =T c 



cir 



dS c 
dT ctf 



dSc 

dm 



(34) 



Using the distribution function Eq. (25) we find, as ex- 



pected, that C c is a measure of the fluctuations of the 
energy and configurational entropy of glassy states. We 
obtain for the configurational heat capacity 



r ( 5F > 

2 eS 



(35) 



6 
a 

CD 



Uj 




r5 c = /c B logiV 



overlap, q" 



FIG. 4: Mean-field theory description of the glassy states 
introduced the concept of an effective potential. Effective 
potential is designed to describe the transition between the 
non-ergodic glassy states and ergodic liquid. 



where (SF) 2 = F 2 - F and {SS c y = S 2 C - S c . Here the 
mean values F and S c are determined by Eq. (33). The 



fluctuations of the configurational entropy and frozen 
state energy are then determined by 



D^P [ip] log 2 P [?/>] 



and 



F 2 = / DipP [tp] F [tp]' 



(36) 



(37) 



respectively. Both quantities can be expressed within the 
replica formalism in terms of a second derivative of F (to) 
with respect to m. For example it follows: 



d 2 1 / —2 

— mF( m ) = --(F 2 -F 2 
am 2 1 V 



(38) 



It is then easy to show that Eq. (35) holds. With the 



introduction of the configurational heat capacity into the 
formalism we have a measure for the deviations of the 
number of metastable states from their mean value. The 
analogy of these results to the usual fluctuation theory 
of thermodynamic variables [79 further suggests that C c 
also determines fluctuations of the effective temperature 
with mean square deviation: 



(39) 



Since C c is extensive, fluctuations of intensive variables, 
like T e ff, or densities, like s c = S c /V, vanish for infinite 
systems. However, they become relevant if one considers 
finite subsystems or small droplets. In the context of 
glasses this aspect was first discussed in Ref. [80 . 
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If the replica theory is marginally stable, i.e. the low- 
est eigenvalue of the fluctuation spectrum beyond mean 
field solution vanishes, it was shown in Ref. [55] that 
T e ff agrees with the result obtained from the generalized 
fluctuation-dissipation theorem in the dynamic descrip- 
tion of mean field glasses [5TJ . Typically, the assumption 
of marginality is appropriate for early times right after 
a rapid quench from high temperatures. In this case 
T c ff > T for T < Ta-, i.e. the distribution function of 
the metastable states is not in equilibrium on the time 
scales where mode coupling theory or the requirement for 
marginal stability applies. Above the Kauzmann tem- 
perature it is however possible to consider the situation 
where T c r — T, i.e. where the distribution of metastable 
states has equilibrated with the external heat bath of the 
system. Since we are interested in the restoration of er- 
godicity for T K < T < T A we use T cS = T. Below the 
Kauzmann temperature the assumption T c g — T cannot 
be made any longer as it implies a negative configura- 



degrees of freedom be governed by the Langevin equation 



tional entropy, inconsistent with the definition Eq. (27 1 



It implies that the glass transition is now inevitable even 
if one could wait for arbitrarily long times. While this is 
likely an idealization of mean field theory, it demonstrates 
that the ageing behavior below and above the Kauzman 
temperature are very distinct. 



B. Dynamical interpretation of the replica 
approach 



Ultimately, glass formation is a dynamic phenomenon. 
It is therefore useful and illustrative to offer a dynamical 
interpretation of the replica formalism of Ref. [16] . Be- 
low we will provide a qualitative set of arguments. A de- 
tailed and quantitative derivation of the formalism, based 
on the dynamic theory of glassy systems by Kurchan and 
Cugliandolo [ST], is given in Ref. [S"o] . 

Consider a system characterized by a field variable 
ip(x,t). Let us assume that certain field configurations 
evolve extremely slowly, and equilibrate at a much later 
times than others. In terms of particle coordinates, a nat- 
ural choice for slow and fast variables would be the fidu- 
cial position Rj (t) and the harmonic vibrations Uj (t), 
according to Xj (i) = (i) + Uj (t). In our field theoretic 
description we assume that slow field configurations are 
characterized by ip (x, t) while we continue to refer to fast 
configurations as <p (x, t) . Let the dynamics of the fast 



,d<p(x,t) _ SS (tp,ip) 



dt 



Sep (x,t) 



+ 7](x,t). 



(40) 



where 



S (ip, il>)=S[<p] + | J dx (ip (x) - iJj [x)f (41) 

contains a weak coupling (g -C 1) between fast and slow 
degrees of freedom, just like in the replica approach dis- 
cussed above. Fast variables ip (x, t) are subject to white 
noise with temperature T 



(r) (x, t) r) (x', t')) = 2TT5 (x - x') S(t- t') . 



(42) 



The equilibrium distribution function for fixed ip is then 
1 



p[(p\tp] = 



exp(-^S (<p, V)) 



(43) 



where Z[tp] = J Dipexp(—/3S[ip,^]). Next we assume 
that the dynamics of the slow degrees of freedom is then 
governed by 



' dt 



5F ty] 



+ r) s (x, t) 



(44) 



5ip (x, t) 

with F [ip] = — T log Z [tp]. Thus, during the entire dy- 
namics of the slow degrees of freedom, the fast modes are 
assumed to be in equilibrium already. It is important to 
assume different noise for the slow variables ip (x,t) and 
for the fast variables ip (x, t). The variables ip (x, t) will, 
in general, not equilibrate with temperature T, but may 
be characterized by an effective temperature T c ff > T, 
i.e. we write: 

( Vs (x, t) Vs {x',t')) = 2T s T cfi S [x - x') 8 (t - if) (45) 
Then, the stationary probability distribution for ip is 

(46) 



p[iP) = ^e X p(-l3 eS F[iP]) = Z ^- 



with m — T/T e Q. The joint distribution of ip and ip can 
now be written as p [ip, ip] — p [ip] p[ip\ip]. It is straightfor- 
ward to show that this distribution gives results identical 
to the replica approach. Above the Kauzmann tempera- 
ture, Tx, equilibration of the slow variables is still possi- 
ble and T approaches T c ff, i.e. m — > 1 



IV. GLASS FORMATION IN UNIFORMLY 
FRUSTRATED SYSTEMS 

To gain detailed insight into the possibility of glass 
formation, we now go back to our model described by 
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and analyze the problem Q with W (cp) — ^j-tp 4 analyt- 
ically within the self consistent screening approximation 
(SCSA) [82] , an approximate treatment that is controlled 
by an expansion in 1/N, where N is the number of com- 
ponents of (p, see Ref. [53] . 

The free energy F(m) of the replicated Hamiltonian is 
given in terms of the regular correlation function G(q) 
and the correlation function F(q) = (f a (q) l p b (~q)) for 
a =/= 6, i.e. between the fields in different replicas. The 
latter corresponds to the Edward- Anderson parameter 
signaling a glassy state. For a system in the universality 
class of the random first order transition we expect that 
below a temperature Ta the system establishes an ex- 
ponentially large number of metastable states and long 
time correlations, characterized by the correlation func- 
tion F(x,x') = T _1 lim^oo (cp(x,t)cp(x',0)). These 
long time correlations occur even though no state with 
actual long range spatial order exists. 

Within the SCSA, the relevant part of the free energy 
F(m) is: 



T 



F(m) = (Trln^T 1 + TrlnX> -1 ), 



m 

which determines the configurational entropy S c 

T dm 



(47) 



. Here, 



■m — l 



Q = (G-F)I + FE 
is the correlation function matrix with I a (, 



(48) 



and 



E a b = 1. The symbol Tr in Eq. (471 includes the trace 



of the replica space and the momentum integration. The 
matrix T> is related to Q via 



v- 



where 



n = (G <g> G - F <g> F)I + (F ® F)E 



(49) 



(50) 



is the generalized polarization matrix. The symbol ® 
denotes a convolution in Fourier space. The replicated 
Schwinger-Dyson equation can be written as 



G^I 



(51) 



where S is the self-energy matrix and Gjj (fc) the Hartree 



correlation function of Eq. (Ill 



Within the SCSA the self-energy has diagonal elements 
£g = 2G(g)Dc and off diagonal elements Hp = 2F®Dp 
in replica space, where Dq and Dp being, respectively 
the diagonal and off-diagonal elements of T>. These equa- 
tions form a closed set of self-consistent equations which 



enable us to solve for G and F, and then determine the 
configurational entropy. To evaluate F(m) we use the 
fact that a matrix of the form A a b = (a — b) S a b + b has 
to — 1 degenerate eigenvalues a t — a — b where i = 1, to — 1 
and and one eigenvalue a m = a + (m — 1) b. This yields 



TrlogA= > / -rloga 



(52) 



This result allows us to evaluate F (to) of Eq. (47). Per- 
forming the derivative of the resulting expression with 
respect to to then yields the configurational entropy 



d d q 



F®F 



(U4T)- 1 + G® G 



(53) 



were, s[x] = —x — ln(l — x). 

The analysis of these coupled equations reveals that 
the self energies Eg and Hp are only weakly momentum 
dependent. Then the impact of Eg is solely to renormal- 
ize the correlation length £. However, the emergence of 
Hp leads to a qualitative change. Dimensional analysis 
reveals that Hp is an inverse length squared, which mo- 
tivates one to introduce the Lindemann length A of the 
glass via A~ 2 = — AqHp. An interpretation of this length 
scale in terms of slow defect motion in a stripe glasses 
was given in Ref. [55] . By inspection of the Dyson equa- 
tion and using the fact that G (k) is strongly peaked at 
the modulation wave number fc one finds for k ~ k that 



F(k )<G(k ). 



(54) 



G (fc) vanishes rapidly away from the peak (as does 
Hp(Jca)G (fco))) and it follows from the same equation, 
[5l] that for large \k — ko\ holds: 



F(k)~-H F (k )G 2 (k). 



(55) 



If a solution for F exists, it is going to be a peaked at 
ko, but smaller and narrower than G. Consequently, if a 
stripe glass occurs, the long time limit of the correlation 
function is not just a slightly rescaled version of the in- 
stantaneous correlation function, but it is multiplied by a 
k dependent function that leads to a qualitatively differ- 
ent behavior for different momenta. Once a glassy state 
is formed, configurations which contribute to the peaks of 
G (k) and F (k), i.e. almost perfect stripe configurations, 
are almost unchanged even after long times. Close to ko, 
F (fc) is solely reduced by some momentum independent 
Debye- Waller factor exp (-D) = F (k ) /G (k ). On the 
other hand, configurations which form the tails of G (fc), 



i.e. defects and imperfections of the stripe pattern, disap- 
pear after a long time since now F (k) <C G (k). The ratio 
of both functions is now strongly momentum dependent. 
F (k) becomes sharper than G (k) because local defects 
got healed in time. The length scale that determines 
the transition between these two regimes is the length A. 
This length can therefore be associated with the allowed 
generalized vibrational motions in a potential minima of 
the complex energy landscape of the system. In analogy 
with structural glasses we therefore call A the Lindemann 
length of the stripe glass. 





ffflW 




In Fig. [6] we show our results for S C (Q) for different 
values of e = £<7 - A self consistent solution of the above 
set of coupled equations corresponds to finding stationary 
points Q* with 



dS c /dQ\ Q=Q , = 0. 



(59) 



time 



Q* = is always a solution with S c (0) = 0. In addition 
there is a locally stable solution for e > sa = 11.770574 
with configurational entropy given as S c (Q\) = 
Vk$K d x 0.0220247 and Q* A = Q* (e A ) = 0.0625528. The 
configurational entropy at the metastable maximum van- 
ishes at e = e K = 13.169625 with Q* K = 0.134271. This 
is the behavior that is generally expected in a system a 
with random first order transition. Upon lowering the 
temperature the ergodic, liquid like state becomes more 
and more correlated, i.e. the correlation length growths. 
Once £ reaches a threshold value £a = |^Zo — 1 - 87Zo 
about twice the typical length scale for field modulations 
Iq = 2n/qo, glassy dynamics sets in as consequence of 
an emergence of exponentially many states. The Kauz- 
man temperature is reached once the liquid state corre- 
lation length grows even further and becomes equal to 
U = If Z ^ 2.1J . 



FIG. 5: Schematic presentation of the defects wandering in 
structural glasses. Here A is an average distance over which 
the defects are manage to wander (Lindemann length) and £ 
is the correlation length (see text). Upper panel illustrates 
the situation when defects are relatively close to each other 
so that they disappear with time and system effectively heals 
itself. Bottom panel depict the opposite situation when wan- 
dering of the defects does not produce healing since they are 
on distances I > A from each other (see Ref. [54] for details.) 

For our subsequent analysis we introduce the dimen- 
sionless quantity Q via 



S F = -Q 



9o 2 



(56) 



i.e. A~ 2 = Qc[q. We are now in a position to perform the 
momentum integration and evaluate the configurational 
entropy as function of Q. It follows that 



S c (Q) = Vk$K d 



Mo (1-2) 7T 



where 



z(Q) = l 



(57) 
(58) 




FIG. 6: Plot of the configurational entropy S C (Q) as a func- 
tion of the dimensionless parameter Q = (l/\qo) 2 , Eq. (4.10) 
for different values of another parameter e — £qo- It is in- 
structive to compare this plot with the schematic plot of Fig. 
4. 

It is also interesting to analyze the penalty for a spa- 
tially varying overlap Q — > Q(x). This will be essen- 
tial for the Landau expansion and the droplet arguments 
that follow below. Performing a gradient expansion of 
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the replica free energy yields the correction 



V. REPLICA LANDAU THEORY 



SF(m) 



1)T / d d x 



(vq)' 



(60) 



where 



«n = 



i^S^ [G{k)G[k+q) - G2[k)] (61) 



The prefactor k^jA\ appears from our definition for Q 



(56 1. For d = 3 and r = 1 the integration can be per- 



formed exactly with the following result 



«2 = ^(fM 5 



(62) 



As was anticipated in Ref. [54], the penalty for a spatial 
deformation becomes small as fco — ^ 0, i.e. in the limit 
where the modulation length is very large. Note that at 
T = Ta or T = Tk, the dimensionless product £fco is 
independent of fco, i.e. varying the characteristic scale 
yields a (T A ) oc k]J 2 . 

Thus, we analyzed an analytically treatable example 
for a random first order transition which enabled us to 
identify the underlying physical mechanism for glassiness 
in a uniformly frustrated system. Once equilibrium corre- 
lations are sufficiently strong and the stripe liquid - stripe 
solid transition is kinematically inaccessible, the system 
is governed by slow glassy dynamics. At this point a 
Lindemann length, A, emerges, which is a length scale 
over which imperfections of the stripe pattern manage to 
wander, Fig. [5| Defects of the perfectly ordered state are 
still abundant as the healing length is finite, i.e. the sys- 
tem enters a random solid state. Local order is however 
established. This may not be identical to the most sta- 
ble crystalline, i.e. long ranged ordered, configuration, 
but rather correspond to configurations that are occur- 
ring with high weight close to the spinodal of the system. 
In the context of supercooled liquids, the natural analog 
of such locally correlated configuration are icosahedral or 
body centered cubic short range order. 

Computer simulations on the model discussed above 
with U3 = 0, suggest that the nucleation barriers of the 
crystalline state are rather low [73] ■ However, the in- 
clusion of asymmetric terms in the interaction potential 
(t*3 =/= 0) will likely strengthen first order transitions to 
a crystalline state and enhance the nucleation barriers. 
This behavior was analyzed in Refs. [56, 64, 65 . 



The analysis of the previous section demonstrated that 
uniformly frustrated systems undergo self generated glass 
transitions with an entropy crisis, of the random first or- 
der transition universality class. These calculations are 
on the mean field level and ignore issues related to dy- 
namical heterogeneity, i.e. variations in the typical time 
scales in spatially nearby regions. The relevance of such 
effect is already evident from the fact that fluctuations 
in the heat capacity become important in finite regions, 
as shown above. To capture these phenomena we need to 
go beyond the mean field level and include droplet for- 
mations. Droplets are not regions with different values 
of the field ip (r) that underlies our analysis. Instead, by 
droplets we mean a spatial variation of the Lindemann 
length, i.e. of the parameter Q used above to character- 
ize the overlap between configurations at distant times. 
Thus, it is useful to consider instead of Q a dynamic field 
Qab (r) that generates the replica partition function, i.e. 



Z (m) = J DQ cxp (-/3H cS [Q] 



(63) 



The mean field solution of this problem is then deter- 
mined by 5H ff [Q] / SQab = with the Ansatz 

Qab = Q* (1 - <U) (64) 
with Q* playing the exact same role as the corresponding 



quantity determined by the condition Eq. ( 59 1 in the 
previous section. 

Formally H c s [Q] can be obtained from the field theory 



(65) 



discussed above by multiplying Z (m) of Eq. ( 32 ) by 
unity 

l =f D f t J] <5(^ b (x,x')-^ a (x)</(x')) = 

x,x' ,ab 
zoo , 

= [ DQ [ Dfi exp S — / ^ Qab ( x ' x ') x 

-L I x ' x ' ab 

x [ Mab (x,x')-^ a (x)^(x')]} 
and switch the order of integration. Then we find 

e -HMQ] = J DfxD m ip.e-^=^ n ^x 

x g- 5Z™b=l f XlX i Qat,(x,x')(p ail (x,x')-<p a (x) ¥ ; b (x')) 



(66) 



We consider the Fourier transform of Q a b (x, x') with re- 
spect to the relative coordinates x — x' and approximate 
it by its value at k — fco . This is consistent with our find- 
ing that the self energies of the mean field theory depend 
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only weakly on momentum. We do however include the 
dependence on the center of gravity coordinate to analyze 
spatially varying overlaps between distant configurations, 
yielding Q ab (x). 

In order to keep our calculation transparent we will 
not analyze the model discussed for uniformly frustrated 
systems, but instead start from a simpler Landau theory 
that is in the same universality class [551 155] : 



[Q] = E 2_, I ^ [Qab] — r 2_j QabQbcQca 

a,b J °0 V c ) 

(67) 



with 

h [Qab] = a \ (VQ afc ) 2 - 



t 2 U- 



:Qab ~ 



II 



Qab + T Qab (68) 



3 

and replica index a,b = 1, . . . , m. is a length scale of 
the order of the first peak in the radial distribution func- 
tion of the liquid and Eq is a typical energy of the problem 
that determines the absolute value of Ts^. In addition 
the problem is determined by the dimensionless variables 
t, u, w and y, which are in principle all temperature de- 
pendent. We assume that the primary T-dependence is 
that of the quadratic term, where t = T ~j^° ■ In what 
follows we further simplify the notation and measure all 
energies in units of E and all length scales in units of 
ao- Formally, Eqs. (67 1 and (68) can be motivated as 



the Taylor expansion of the free energy F [m] with re- 
spect to S a f,. In practice, the explicit Taylor expansions 
of more microscopic models do not yield quantitatively 
reliable results. For example, the Taylor expansion of 



Eq. (57) poorly reproduces the behavior in the glassy 



regime. Still, the simple Landau expansion of Eqs. (67) 



and ( 68 1 allows us to gain key insights that enable us to 
investigate the more general case. In particular, we will 
demonstrate that the replica instanton theory can easily 
be generalized to the more complex stripe glass model. 
The mean field analysis of the model Hamiltonian Eqs. 



( 67 1 and ( 68 ) is straightforward. Inserting the Ansatz Eq. 



(64) into H c ff [Q] and minimizing w.r.t Q yields Q* = 
or 



Q* (t) 



4ty 



2/y 



(69) 



Nontrivial solutions exits for t < tA — ^ with Q* A = 
Q* (^a) = f^j which determines the mode coupling tem- 
perature T A = t A E + T - Inserting Q* of Eq. (|69| into 
H c s [Q] yields at the saddle point F (m) = H c s [Q] for 
the replica free energy, which yields: 



F(m) 
V(m- 1) 



t 



*2 



+ u (m — 1) 



1*3 



3' 



(70) 



The configurational entropy, as determined by S c 



k-Jr-F {m)\ , is given by 

T am \ / 1 m— >1 & J 



Sc(Q) 



v ft 



T \2 



,Q 2 



w 
— ( 

3 



(71) 



The stationary points of S c (Q) are given by Q* of Eq 69 

2w' 
9y 



Inserting Q* yields that S c vanishes at tx = with 



Q*k = 1^- Close to tfc it follows that 



V- 



tx (t — tx 
~Tk~ 



oc V 



T-T K 



as expected. At t A one finds T A S C (T A ) = V ig 2y 3 
the configurational heat capacity follows 



C r . = 



2V 



\Q* 2 - ^pQ* 3 + |<3* 4 



(72) 
For 

(73) 



It holds at C c (Trt) = 



Vu 2 
T K 3 



S c 



■ill 



C c {t 



t-t K 



K 



t K 



(74) 



Thus, we see that the main findings of the model calcu- 
lations are reproduced by the simple Landau expansion, 



Eqs. (67) and (68). They will now be used as starting 



point for our analysis of dynamical heterogeneity in form 
of a replica instanton theory. 



VI. REPLICA INSTANTONS AND ENTROPIC 
DROPLETS 

At the mean field level a glass at T < Ta is frozen 
in one of many metastable states. Q* then characterizes 
the overlap between configurations at distant times. The 
free energy of such a frozen state is higher by TS C com- 
pared to the ergodic liquid state that is characterized by 
Q* = 0. Thus, for T K < T < T A the mean field glass 
is locally stable. Local stability also follows from the 
fact that the lowest eigenvalue of the fluctuation matrix 
6 2 H/SQ ab (x) 5Q cd (x') is positive for T K < T < T A and 
vanishes at T = Ta, see Ref. [58] . The Q-dependence 
of S c (Q) shown in Fig. [6] suggests that the decay modes 
for the frozen state are droplet excitations, similar to the 
nucleation of an unstable phase close to a first order tran- 
sition. This situation was analyzed in Ref. [55] [53]. In 
agreement with the RFOT theory [2D], the driving force 
for nucleation is the configurational entropy, leading to 
the notion of entropic droplets. The formal approach to 
analyze entropic droplets is performed in terms of the 
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effective potential approach of Refs. [HO [EH [SS]. We 
used this technique to formulate a replica instanton and 
barrier fluctuation theory in Refs. [23 In what fol- 
lows we use a slightly simpler approach that yields es- 
sentially the same results, but is physically significantly 
more transparent. 

As can be motivated by the more general effective po- 
tential approach of Refs. [19, 85, 86 , instanton solutions 
for entropic droplets can be determined from 



SH [Q] 

SQab (x) 



= 0, 



(75) 



where we allow for spatial variations of the overlap 
Qab ( r ) = Q ( r ) (8 a b ~ !)• This yields the following non- 
linear equation for Q (r) : 



V 2 Q(r) 



ds c (Q (r)) 
dQ (r) : 



(76) 



with configurational entropy density s c (Q) = S c (Q) /V . 
In case of the Landau expansion we use S c (Q) of Eq. 



(711, while for the stripe glass approach we start from 
Eq. (57). We will first analyze the simpler case of the 



Landau theory. Then Eq. ( 76 ) admits an exact solution 
in the thin wall limit R ^> £: 



Q(x) 



th 



ZQ 



th 



(77) 

where the integration constant zq is a function of i, w and 
y. R is the droplet radius and £ is the interface width 
given by 



4a n 



v/3y(2Q* -Q* K y-m K + At 



(78) 



Inserting the solution Eq. ( |77j ) into the expression into 
H [Q] we calculate the value of the mean barrier. The 
latter is determined by optimizing the energy gain due 
to creation of a droplet and energy loss due to the sur- 
face formation. As a result for the mean barrier we find 
(reintroducing the energy scale Eq and length scale do) 



^ 327roo 2 
9^ 3 ' 



The droplet radius 
R = 



64ag 



3y 2 q* 3 (q* K -q* (t))£ 3 



(79) 



(80) 



is determined from the balance between the interface ten- 
sion and the entropic driving force for nucleation. Fur- 
thermore, q* K = q* (t — tx) is the order parameter at the 



Kauzmann temperature. When temperature approaches 

the T K the radius of the droplet as well as the mean 

barrier diverge. One finds lim Ft oc (t — tK)~ 2 and 

t— Hk 

lim R oc (t — t^-) -1 . Since the droplet interface £ re- 

t— >tK 

mains finite as t — > fx, the thin wall approximation is 
well justified close to the Kauzman temperature. On the 
other hand, R and £ become comparable for tempera- 
tures close to Ta and the thin wall approximation breaks 
down. Combining R oc (t — tx) and s c oc (t — tx), we 
obtain v = 1 for the exponent that relates the droplet 
size R and the configurational entropy density: R oc s~ v . 
Close to t K follows <t (t) ~a K (l-^ (^tjr)) with 
§j2 . The critical droplet nucleation radius is 



a K = 



4ap w 



27v / 15 V 



R = 2a/s c yielding the barrier 



Ft 



167T(7 J 
3 S 2 



(81) 



for the nucleation of entropic droplets. This leads to a 
mean relaxation time of 



r exp 



Ft 



(82) 



It was pointed out in Ref. [50] that wetting effects 
of the interface alter the a relationship between droplet 
size and entropic driving force to R oc s~ u with expo- 
nent v. A renormalization of the droplet interface due to 
wetting of intermediate states on the droplet surface was 
shown to yield v = 2/d [5D], leading to Ft oc Ts~ x and 
correspondingly to a Vogel-Fulcher law 



T exp 



DT K 
T-T, 



K 



(83) 



for the mean relaxation time. 

It is straightforward to analyze the more complex con- 
figurational entropy of the stripe glass problem. In this 
case the nonlinear instanton equation cannot be solved 
analytically. However we can make progress by perform- 
ing a variational calculation for the droplet. The mean 
barrier is 



Ft = H sS [Q (x)]-5 c (Q*) s 



(84) 



where Q (x) is a localized instanton solution which differs 
from Q* in a finite region. We make the trial ansatz: 



Q(x) 



Q*\JS\^R R<\ X \<R 



x| < R 
|x| < iU 
x| > R + l 



(85) 
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and insert it into the above expression for H c s [Q (x)]. 
We find 



Ft (R) = AiraR 2 



-in 



s c (Q*)R 3 



with surface tension 



1*2 



21 



IT 



8c (Q) 



dQ 



s c (Q* 



(86) 



(87) 



Minimizing with respect to the droplet wall thickness 
yields 



I 2 



1 



a\Q 



*3 



2 s c (Q)dQ-s c (Q*)Q* 
and correspondingly for the surface tension 

a = V2a Q* 1/2 (j^ s c (Q) dQ - s c (Q*) Q* 



1/2 



(89) 



Using the our result Eq. (57) for the configurational en- 



tropy as a function of Q for the stripe glass problem in 
d = 3 and for r = 1 (i.e. with the long range Coulomb in- 



teraction V(x)) and Eq. (62) for the gradient coefficient 
ao, yields at Tk ■ 

°{Tk) = C a T K kl (90) 

where C a ~ 3.45 x 10~ 2 . Thus we see that the surface 
tension of entropic droplets vanishes in the limit ko — > 0. 
For the wall thickness we obtain 1(Tk) — k^ 1 . 

Finally, we use our results for the surface tension to 



analyze the variation of the fragility D in Eq. (83) as 
function of fco . One finds 



D 



3<r 2 



(91) 



T=T K 



which was further investigated in a numerical analysis of 
the stripe glass problem in Ref. 90 . This analysis led 
to D(ko — > 0) —> 0. If we make the following estimate 



dS c 
~dT 



T=T R 



S C (T A 



K 



(92) 



and neglect the dependence of Ta — Tk on ko, we find 
that D oc fco in qualitative agreement with the results of 
Ref. M. 



A. barrier fluctuations 

Numerous experiments on supercooled liquids are not 
only sensitive to the mean barrier, F$, but are able to 



measure the entire (broad) excitation spectrum in glasses 
|87j . Most notably, the broad peaks in the imaginary 
part of the dielectric function e" (u) are most naturally 
understood in terms of a distribution g (r) of relaxation 
times, such that 



e" (w) cx / drg (r) 



2 ' 



(93) 



Similarly dynamical heterogeneity with spatially fluctu- 
ating relaxation times yields non-exponential (frequently 
stretched exponential) relaxation of the correlation func- 
tion 



•(*) 



drgir) £ - t/T , 



(94) 



Other effects that are most likely caused by a distribu- 
tion of relaxation rates include the break down of the 
Stokes-Einstein relation D = ^ bT t between the diffusion 
coefficient D of a particle of size L and the viscosity r\ 
[88l [89] . These experiments call for a more detailed anal- 
ysis of the fluctuations 



5F$ 2 = Ft 2 - Ft 



(95) 



of the activation barriers and, more generally, of the dis- 
tribution function p (F*) of barriers. The latter yields 
the distribution function of the relaxation times 



g(r)=p{Fi) 



dF* 



(96) 



through r (F*) = r exp ( ^j^j ■ For example, in case of 
a Gaussian distribution of barriers one obtains a broad, 
log-normal distribution of relaxation rates: 



•(r) 



1 



: exp 



log 2 (t/t) 
2A 



with 



A = 5Ft 2 / (hsT) 



(97) 



(98) 



and t from Eq. (82 1. While the distribution, Eq. <\97l, 



does not yield a stretched exponential form for the cor- 
relation function, it can often be approximated by 



(99) 



with (3 = (1 + A) ' . Furthermore, the study of higher 
order moments of p (F*) is important to determine 
whether the distribution is indeed Gaussian or more com- 
plicated. 

In Ref. |59j we used the replica formalism discussed 
here as well as the more elaborate replica method of 
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formalism Refs. [TH1 [EH US] to determine higher order 
moments of the barrier distribution function. Using the 



"thin wall" approximation for Q(x) given by Eq. (77) 
we obtain for the second moment 



5F& =A(R i + P R 2 



(100) 



where R is the radius of the droplet where the explicit 
expressions for the coefficient A and length p are given 
in Ref. [SS]. It is noteworthy that there is a surface 
contribution to the moment of the barrier fluctuations 
that is a consequence of correlations between droplet and 
homogeneous background with overlap Q* . In Ref. [SH] 
it was also shown that the barrier distribution function 
p (-F^) is Gaussian, at least if one considers the first six 
moments. The skewness of the actual barrier distribution 
measured in structural glasses is, in our view, an effect 
due to the interaction of spatially overlapping instantons. 



VII. SUMMARY 

To summarize, we have discussed several models of 
dassy systems where the randomness is self generated 



rather than induced by strong external factors. In par- 
ticular, we applied the replica formalism developed for 
the spin glass systems to study the glass transition in the 
many-body systems at the presence of an arbitrary weak 
disorder. We employed the Landau theory to analyze the 
mean field glass transition using the saddle point approx- 
imation. We have also considered the energy fluctuations 
around the saddle point and evaluated the barrier height 
distribution. 
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